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Abstract Several plant pathogenic Parastagonospora species have been identified infecting wheat and other 
cereals over the past 50 years. As new lineages were discovered, naming conventions grew unwieldy and the 
relationships with previously recognized species remained unclear. We used genome sequencing to clarify relation- 
ships among these species and provided new names for most of these species. Six of the nine described Parasta- 
gonospora species were recovered from wheat, with five of these species coming from Iran. Genome sequences 
revealed that three strains thought to be hybrids between P. nodorum and P. pseudonodorum were not actually 


hybrids, but rather represented rare gene introgressions between those species. Our data are consistent with the 
hypothesis that P. nodorum originated as a pathogen of wild grasses in the Fertile Crescent, then emerged as a 
wheat pathogen via host-tracking during the domestication of wheat in the same region. The discovery of a diverse 
array of Parastagonospora species infecting wheat in Iran suggests that new wheat pathogens could emerge from 


this region in the future. 
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INTRODUCTION 


Parastagonospora nodorum is an important wheat pathogen 
with a global distribution. It has been subjected to intensive 
population genetic analyses using a variety of genetic markers, 
including RFLPs in nuclear and mitochondrial genomes (Keller 
et al. 1997, Sommerhalder et al. 2007), microsatellites (SSRs; 
Stukenbrock et al. 2006), sequences of genes encoding impor- 
tant traits like virulence (Stukenbrock et al. 2007b, Ghaderi et al. 
2020) and fungicide sensitivity (Pereira et al. 2017) as well as 
entire genome sequences (Richards et al. 2019, Pereira et al. 
2020a, c). These studies revealed that E nodorum populations 
are characterized by regular recombination, high levels of gene 
flow, high effective population sizes, and varying frequencies of 
effector genes encoding host-selective toxins (McDonald et al. 
2013). Analyses of several field populations from Iran indicated 
that P. nodorum most likely originated in the Fertile Crescent 
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(McDonald et al. 2012, Ghaderi et al. 2020), the same region 
where wheat was domesticated (Salamini et al. 2002), providing 
support for the hypothesis that this pathogen emerged through 
host-tracking (Stukenbrock et al. 2008). Analyses of DNA 
sequences revealed that several different Parastagonospora 
lineages can be found on wheat in Iran, including one spe- 
cies that was also found infecting wheat in other geographical 
regions, as well as new lineages that had not previously been 
found on wheat (McDonald et al. 2012). Some lineages were 
also suspected to be of hybrid origin (McDonald et al. 2012), 
which suggested there may be significant uncertainty associ- 
ated with resolving phylogenetic relationships among these 
lineages based on a small number of gene sequences. Here, we 
used genome sequences to identify thousands of orthologous 
gene sequences to build a phylogenomic tree among these 
wheat-infecting species and place them into a larger context 
that includes Parastagonospora species found on wild grasses 
in the Fertile Crescent and other continents (Appendix). 


The different pathogens in the Parastagonospora complex (pre- 
viously Phaeosphaeria) causing cereal diseases were originally 
defined based on spore morphology, host specialization and 
the formation of sexual ascomata (Shaw 1957a, b). Isolates 
that were most aggressive on wheat and showed heterothallic 
mating type behaviour formed the group now called P. nodorum. 
Isolates collected from oats or other hosts were initially named 
Leptosphaeria avenaria, with the species name reflecting the 
host preference. Subsequent DNA-based analyses supported 
placing this group into the Phaeosphaeriaceae and it was re- 
named Phaeosphaeria avenaria. A third group of isolates were 
morphologically similar to P. avenaria, but were not able to infect 
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oats, instead being weakly aggressive on wheat. This group 
was named Phaeosphaeria avenaria f. sp. tritici (Pat) (Ueng et 
al. 1995). Further analyses that provided additional evidence 
for host specialization (Martin & Cooke 1979, Osbourn et al. 
1986, Cunfer & Ueng 1999) and genetic differences between 
host specialized forms (Ueng & Chen 1994, Ueng et al. 1998, 
Malkus et al. 2005) led to splitting Pat into three groups called 
Pat-1, Pat-2 and Pat-3. These early phylogenetic studies, which 
often relied on analyzing a single genetic locus and a small 
number of isolates, gave inconsistent results. 


Amore recent study (McDonald et al. 2012) that included more 
than 350 globally distributed Parastagonospora strains and 
sequences of five different genes (both mating type loci, ITS, 
B-tubulin and B-xylosidase) led to the discovery of new Pat line- 
ages infecting other grass species in Iran and North America. 
This dataset also identified 14 strains in Canada and Iran that 
appeared to be hybrids between Pat 7 and P nodorum, sug- 
gesting that species boundaries may be permeable. The Iranian 
wheat fields were colonized mainly by P. nodorum and Pat-1, but 
also by two other Parastagonospora species identified earlier 
(McDonald et al. 2012) and a new species discovered in this 
analysis. Arecent study used ITS sequences to identify a clade 
of Parastagonospora strains infecting native Stipa pulchra grass 
in California (Spear & Mordecai 2018). The affiliation of this 
clade to the other Parastagonospora species found on grasses 
was not known, but we hypothesized that this species may 
have emerged through the mechanism of ‘spill-back’, where a 
pathogen specialized to infect a host — often a domesticated 
host — in one location is introduced into a new location and 
becomes adapted to a new host — often a native, wild species 
(Kelly et al. 2009). An example of spill-back in plant pathology 
is the adaptation of Rhynchosporium graminicola to infect wild 
barley (Hordeum spontaneum) in northern Africa thousands of 
years after its emergence as a pathogen of cultivated barley 
(Hordeum vulgare) in northern Europe (Kiros-Meles et al. 2011). 


We utilized a genome-scale dataset and a wide array of Para- 
stagonospora species to explore the phylogenetic relationships 
among these species and to better define the species bounda- 
ries. Based on our findings, we propose names for seven 
new Parastagonospora species. We also sought to determine 
whether the putative Pat-1/P. nodorum hybrid strains identified 
based on analyzing five genes were truly hybrids. Finally, we 
explored the hypothesis that the Parastagonospora strains 
found on wild Stipa could have originated via spill-back from 
P. nodorum strains infecting wheat in California. 


MATERIALS AND METHODS 


Isolates 


Twenty-two of the analyzed Parastagonospora strains were 
isolated from wheat, including four strains originating from 
global collections of P. nodorum that were already analyzed 
for both phenotypes and entire genome sequences (Pereira 
et al. 2020a, b, c). Seven strains came from other grass hosts 
and had previously been analyzed using a limited number of 
conserved gene sequences (McDonald et al. 2012), but had 
not previously been genome sequenced. We included repre- 
sentatives of most of the known Phaeosphaeria avenaria f. sp. 
tritici (Pat) clades that were identified by plant pathologists over 
the last three decades. Parastagonospora stipae isolates were 
obtained from wild Stipa pulchra in the Jasper Ridge Biologi- 
cal Preserve, San Mateo County, California during a survey 
of naturally occurring foliar fungal pathogens in the summer 
of 2015 (Spear & Mordecai 2018). A summary of the isolates 
included in these analyses is shown in Table 1. 


For strains included in morphological analyses, single conidial 
colonies were established by excising 2 x 2 mm tissue pieces 


from the leading margin of disease, surface sterilizing them in 
sequential baths of alcohol and bleach, and plating it onto 2 % 
malt extract agar (MEA) with chloramphenicol. These isolates 
were then sub-cultured to obtain axenic cultures. Cultures 
were maintained in Petri dishes sealed with parafilm at room 
temperature in the laboratory. Colonies were sub-cultured on 
2 % potato-dextrose agar (PDA), oatmeal agar (OA), MEA, 
autoclaved Triticum and Hordeum seed and leaves on synthetic 
nutrient-poor agar (SNA) (Crous et al. 2019), and incubated 
at 25 °C under continuous near-ultraviolet light to promote 
sporulation. Reference strains and specimens of the studied 
fungi are maintained in the CBS culture collection (CBS) of 
the Westerdijk Fungal Biodiversity Institute (WI), Utrecht, the 
Netherlands. 


Culturing, DNA extraction, and whole-genome sequencing 


All isolates were initially established from single-spore isola- 
tions. Four-day-old colonies growing on PDA media were har- 
vested for mycelial fragments and transferred to 50 mL Potato 
Dextrose Broth (PDB) media. The transferred fragments were 
cultured 4—6 d at 24 °C at 120 rpm. Sterile cheesecloth was 
used for filtering and the fungal material was then lyophilized 
at room temperature for 72 h. DNAextractions were performed 
from dried material with the DNeasy Plant Mini Kit (Qiagen) 
following the manufacturer’s protocol. Genomic DNA was se- 
quenced on the Illumina HiSeq 2500 platform with a paired-end 
100 bp cycle protocol. Library preparation and sequencing was 
performed by the Functional Genomics Center Zurich (FGCZ). 
Raw sequence reads were deposited in the NCBI Sequence 
Read Archive (SRA) under the BioProject PRJNA606320 
(Pereira et al. 2020a, c). Sequencing reads of additional ge- 
nomes were retrieved from BioProject accessions SRP155908, 
SRP159197 (Syme et al. 2018, Richards et al. 2019). 


Genome assembly and gene annotation 


Draft genome assemblies were produced for each isolate us- 
ing SPAdes v. 3.14.1 (Bankevich et al. 2012) with the --careful 
option and a pre-determined kmer range of 21,33,55,66,99,127 
(the latter kmer was omitted if the read length was 100 bp). As- 
semblies were quality checked using Quast v. 5.0.2 (Gurevich et 
al. 2013). To train accurate gene predictions, RNAseq data from 
the P. nodorum reference genome strain SN15 was retrieved 
from the NCBI Sequence Read Archive (Jones et al. 2019). 
The datasets SRR11785359, SRR11785360, SRR11785362 
matching the wild-type genotype of SN15 were used. All 
transcriptomic reads were aligned to the reference genome 
of P. nodorum SN15 (Hane et al. 2007) using STAR v. 2.7.5a 
(Dobin et al. 2013). bam2hints included in Augustus v. 3.3.3 
was used to retrieve intron splice site hints and filtered for a 
minimum read support of 10 (Stanke & Morgenstern 2005). 
Intron splice site hints were used for a comprehensive gene 
model training using Braker v. 2.1.5 including GeneMark-ET 
and Augustus v. 3.3.3 (Stanke & Morgenstern 2005, Hoff et al. 
2019). After successful gene model training, Augustus v. 3.3.3 
was used to annotate all assembled genomes with ‘--alterna- 
tives-from-evidence=false’ and ‘--UTR=off’. Predicted protein 
sequences were retrieved from each genome for downstream 
analyses. 


Orthology analyses of maximum-likelihood tree recon- 
struction 


Orthology relationships were established among all assembled 
genomes as well as genomes from representative ascomycetes 
retrieved from Fungal Ensembl Genome (release 47) using a 
script provided by G. Leonard (https://github.com/guyleonard/ 
get_jgi_genomes). Sets of protein sequences were used for 
pairwise BLAST analyses and ortholog reconstruction using 
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Orthofinder v. 2.4.0 (Emms & Kelly 2019). The default infla- 
tion index of 1.5 was used and set the sequence alignment 
procedure to ‘msa’, which uses mafft with --maxiterate 1000 
to generate sequence alignments of each orthogroup (Katoh & 
Standley 2013). Aconcatenated alignment was retrieved of all 
single-copy orthogroups (n = 2425) with a minimum of 98.0 % of 
species having single-copy genes in any orthogroup. Maximum 
likelihood tree searches were performed with a randomized 
parsimony starting tree, a general time reversible (GTR) protein 
model with unequal rates and unequal base frequencies us- 
ing raxml-ng 0.9.0 (Kozlov et al. 2019). Felsenstein bootstrap 
replicates (n = 100) were performed and the final tree was 
produced using ggtree (Yu 2020). 


Sequencing of fungal barcoding genes, cleaning, 
concatenation and phylogeny 


Primers for the AFTOL (Assembling the Fungal Tree of Life) 
genes (ITS, LSU, RPB1, RPB2, TEF1) were designed using 
the software Primer3 (Untergasser et al. 2012) based on the 
reference genome of the SN15 isolate (Hane et al. 2007). The 
genes were amplified in P. stipae isolates with conditions for 
PCR amplification as follows: 96 °C for 2 min, 35 cycles at 
96 °C for 30 s, 56 °C for 30 s, 72 °C for 1 min, and a final step 
at 72 °C for 5 min. Amplified products were purified using manu- 
facturer protocols for illustra Sephadex G-50 fine DNA Grade 
Column (GE Healthcare, Pittsburgh, USA) and sequenced in 
an ABI 3130xI Genetic Analyzer (Life Technologies, Applied 
Biosystems). Raw Sanger sequence reads of each primer were 
checked for quality and assembled into final gene sequences 
using Geneious v. 9.1.8 (Biomatters, Auckland, New Zealand). 
Gene sequences for ITS, LSU, RPB2 and TEF1 were available 
for P. phragmitis and P. novozelandica (Marin-Felix et al. 2019). 
Nucleotide sequences for the remaining species were extracted 
from the genome sequences available. Final datasets were 
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obtained after a concatenation step using the MAFFT online 
tool v. 7.0 (Katoh et al. 2019). The concatenated dataset files 
were used to reconstruct maximum likelihood trees also in 
RaxML v. 8.0 (Stamatakis 2014). 


Morphology 


Representative strains were morphotyped for a series of mor- 
phological characters. Slide preparations were mounted in 
Shear’s mounting fluid from colonies sporulating on MEA, 
PDA, OA or leaf tissue on SNA. Observations were made 
with a Nikon SMZ25 dissection-microscope, and with a Zeiss 
Axio Imager 2 light microscope using differential interference 
contrast (DIC) illumination and images recorded on a Nikon 
DS-Ri2 camera with associated software. Colony characters 
and pigment production were noted after 2—4 wk of growth 
on MEA, PDA and OA (Crous et al. 2019) incubated at 25 °C. 
Colony colours (surface and reverse) were scored using the 
colour charts of Rayner (1970). 


RESULTS 


The genome-scale tree resolved nine distinct Parastagono- 
spora species (Fig. 1). Parastagonospora nodorum isolates 
from wheat fields on four continents (Africa, North America, 
Europe, Australia) formed a single well-resolved clade that 
was most closely related to P. pseudonodorum (previously 
named Pat-1), P. bromicola (previously Pat-5), P. dactylidigena 
(previously Pat-4) and P. avenae (previously Phaeosphaeria 
avenaria). Considerable differences were found between the 
two P. avenae strains, which may reflect additional cryptic 
species in this clade. More distant from P. nodorum were 
P. golestanensis (previously Pat-6), P. zildae (previously P-1), 
P. jasniorum (previously P-2) and P. arcana. Five of these spe- 
cies (P. arcana, P. jasniorum, P. nodorum, P. pseudonodorum, 
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Fig. 1 Maximum likelihood phylogenomic tree generated from a concatenated alignment of sets of orthologous protein sequences retrieved from Parastago- 
nospora draft genomes and representative ascomycetes. A total of 2425 single-copy orthologs with = 98.0 % of the species included were retained. The tree 
was estimated based on the general time reversible (GTR) protein model and branch support corresponds to Felsenstein bootstrap values (n = 100). Bootstrap 
values below 98 % were omitted. The root was defined as the node connecting Tuber melanosporum. Genome sequences outside of the Parastagonospora 
genus were retrieved from Ensembl Fungi (https://fungi.ensembl.org/index.html). 
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Fig. 2 Maximum likelihood phylogenetic trees. a. Phylogenetic tree generated from a concatenated alignment of ITS, LSU, RPB1, RPB2 and TEF1a nucleotide 
sequences that include Parastagonospora stipae; b. phylogenetic tree generated from a concatenated alignment of ITS, LSU, RPB2 and TEF1 sequences 


that includes Parastagonospora phragmitis and P. novozelandica. 


P. zildae) were found infecting wheat in Iran and P. avenae was 
found infecting wheat in North America, indicating that at least 
six species in the Parastagonospora clade can infect wheat. 


Because we did not obtain genome sequences for any of the 
P. stipae isolates, we constructed a new phylogenetic tree using 
concatenated ITS, LSU, RPB1, RPB2 and TEF 1a sequences for 
all of the Parastagonospora species (Fig. 2a). Many of the spe- 
cies relationships indicated in the genome-scale phylogenetic 
tree were also found in the tree based on these five AFTOL 
genes. This tree indicated that P. stipae was distinct from the 
other species, but was most closely related to P. avenae. 


Parastagonospora phragmitis and P. novozelandica are two 
recently described Parastagonospora species infecting wild 
grasses in Australia and New Zealand, respectively (Marin-Felix 
et al. 2019). Their relationship to the nine Parastagonospora 
species described here was determined based on concaten- 
ated ITS, LSU, RPB2, TEF1a sequences (Fig. 2b). Many 
of the species relationships indicated in the genome-scale 
phylogenetic tree were also found in the tree based on these 
four AFTOL genes. This tree indicated that P. phragmitis and 
P. novozelandica were separate species that were most closely 


related to P. golestanensis, while P. stipae was most closely 
related to P nodorum and a P. avenae strain isolated from 
wheat in New York. 


Taxonomy 


Based on morphology and a multi-locus phylogeny of the iso- 
lates studied, nine Parastagonospora species are recognized, 
of which seven are newly described. In total, 26 species of 
Parastagonospora are presently known from DNA sequence 
data (Appendix). 


Parastagonospora arcana B.A. McDonald, P.C. Brunner, Croll, 
D. Pereira & Crous, sp. nov. — MycoBank MB838010; Fig. 3 


Etymology. arcanus = secret or mystery, reflecting the fact that this spe- 
cies was previously assumed to be P. nodorum. 


Typus. Iran, Golestan Province, on infected glume of Triticum aestivum, 
2005, M. Razavi (holotype CBS H-24479, cultures ex-type CPC 36218 = 
IR_G2.1A = CBS 146965). 


Ascocarps pseudothecial, solitary, erumpent, globose, brown to 
dark brown, 200—250 um diam; wall of 3—6 layers of brown 


Fig. 3 Parastagonospora arcana (CPC 36221). a. Ascomata developing on PDA; b—d. asci, ascospores and pseudoparaphyses. — Scale bars: a = 250 um, 


all others = 10 um. 
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textura angularis; ostiole central, not to slightly papillate. Asci 
70-95 x 8—10 um, bitunicate, subcylindrical, straight to slightly 
curved, short stipitate, 8-spored, with well-defined apical cham- 
ber, 1-1.5 um diam. Pseudoparaphyses filiform, hyaline, sep- 
tate, hyphae-like, unbranched, rarely anastomosing, 2—2.5 um 
diam. Ascospores bi- to triseriate in asci, fusoid, guttulate, 
pale brown, smooth, widest in penultimate cell, becoming 
slightly constricted at septa with age, (24—)25—28(—32) x (4—) 
5(—5.5) um, 5-septate. 

Cultural characteristics — Colonies covering dish after 2 wk 
at 25 °C in the dark, smooth, with even margin, and moder- 
ate to fluffy aerial mycelium. On MEA surface pale olivaceous 
grey, reverse isabelline; on PDA surface pale olivaceous grey, 
reverse isabelline; on OA surface pale olivaceous grey. 


Notes — Parastagonospora arcana is sister to P. zildae 
(98.4 % identity based on ITS) and confirmed as distinct based 
on the multigene dataset. Both species only produced a sexual 
morph in culture. Morphologically, the two species are very 
similar, and are best distinguished based on their DNA data. 


Parastagonospora avenae (A.B. Frank) Quaedvlieg et al., 
Stud. Mycol. 75: 362. 2013 — Fig. 4 


Basionym. Septoria avenae A.B. Frank, Ber. Deutsch. Bot. Ges. 13: 64. 
1895. 

Synonyms. Stagonospora avenae (A.B. Frank) Bissett (as ‘avena’), Fungi 
Canadenses, Ottawa 239: 1. 1982. 

Leptosphaeria avenaria G.F. Weber, Phytopathology 12: 449. 1922. 

Phaeosphaeria avenaria (G.F. Weber) O.E. Erikss., Ark. Bot., Ser. 2 6: 
408. 1967. 

Pleospora tritici Garov., Arch. Triennale Lab. Bot. Crittog. 1: 123. 1874. 


Spermatogonia solitary to aggregated, globose, brown, 60—90 
um diam with central ostiole; wall of 3—6 layers of brown textura 
angularis. Soermatophores reduced to spermatogenous cells. 


Spermatogenous cells hyaline, smooth, ampulliform, 7—10 x 
3—5 um, with indistinct apical percurrent proliferation. Soermatia 
solitary, hyaline, smooth, guttulate, subcylindrical with obtuse 
apex and truncate base, aseptate, 4—6 x 2 um. 

Culture characteristics — Colonies flat, spreading, reaching 
60 mm diam after 2 wk, with moderate aerial mycelium, and 
even, smooth margins. On MEA, PDA and OA surface and 
reverse ochreous. 


Material examined. Austra, Mt Barker, on Avena sativa, 2009, K. Clarke, 
specimen CBS H-2439, culture CPC 36201 = avenaria_BARKER = CBS 
146876. 


Notes — Bissett (1982) noted that conidia in vivo can be (1—) 
3(-7)-septate, 17—46 x 2.6—4.4 (av. 33 x 3.5) um in size. Un- 
fortunately, no conidiomata were observed in culture as the 
present strain only formed spermatogonia. 


Parastagonospora bromicola B.A. McDonald, P.C. Brunner, 
Croll, D. Pereira & Crous, sp. nov. — MycoBank MB838011; 
Fig. 5 


Etymology. Name reflects the host genus it was isolated from, Bromus. 


Typus. USA, North Dakota, on Bromus inermis, 1983, J. Krupinsky, 
USDA-ARS (holotype CBS H-24387, cultures ex-type CPC 36214 = Pat-5 
83.6011.2 = CBS 146870). 


Conidiomata solitary, pycnidial, dark brown to black, globose, 
150-200 um diam, with central ostiole; wall with 3—6 layers 
of pale brown textura angularis. Conidiophores reduced to 
conidiogenous cells lining the inner cavity. Conidiogenous cells 
hyaline, smooth, ampulliform to subcylindrical, 4—6 x 4—5 um, 
proliferating percurrently at apex. Conidia solitary, hyaline, 
smooth, guttulate, subcylindrical, straight to slightly curved, 
apex subobtuse, base truncate, 2 um diam, 1(—3)-septate, (12—) 
14-16 (-18) x (2—)2.5-3 um. 


Fig.4 Parastagonospora avenae (CPC 36201). a. Conidiomata developing on PDA; b. conidiomal wall with conidiogenous cells; c. close-up of conidiogenous 
cells giving rise to conidia; d. conidia. — Scale bars: a = 100 um, b = 50 um, c—d = 10 um. 


Fig. 5 Parastagonospora bromicola (CPC 36214). a. Conidiomata developing on PDA; b. conidioma with oozing conidia; c. conidiogenous cells giving rise 


to conidia; d. conidia. — Scale bars: a = 200 um, all others = 10 um. 
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Culture characteristics — Colonies flat, spreading, cover- 
ing dish after 2 wk, with moderate aerial mycelium, and even, 
lobate margins. On MEA surface luteous, reverse ochreous; on 
PDAsurface smoke grey, outer region olivaceous grey, reverse 
olivaceous grey; on OA surface saffron. 


Notes — Parastagonospora bromicola occurs on Bromus, 
and has smaller conidia than are commonly associated with 
taxa in the P. nodorum complex. 


Parastagonospora dactylidigena B.A. McDonald, P.C. Brunner, 
Croll, D. Pereira & Crous, sp. nov. — MycoBank MB83801 2; 
Fig. 6 


Etymology. Name reflects the host genus it was isolated from, Dactylis. 


Typus. IRAN, Golestan Province, on Dactylis glomerata, 2011, M. Razavi 
(holotype CBS H-24386, cultures ex-type CPC 36213 = Pat-4 IR_2_5.2 = CBS 
146869). 


Conidiomata solitary, pycnidial, brown, globose, 250—350 um 
diam, with central ostiole; wall with 3—6 layers of pale brown 
textura angularis. Conidiophores reduced to conidiogenous 
cells lining the inner cavity. Conidiogenous cells hyaline, 
smooth, ampulliform to subcylindrical, 5—7 x 4—5 um, prolifera- 
ting percurrently at apex. Conidia solitary, hyaline, smooth, 
guttulate, subcylindrical, straight to slightly curved, apex sub- 
obtuse, base truncate, 2.5—3.5 um diam, 3(—6)-septate, (25-) 
30—37(—42) x 4(—5) um. 

Culture characteristics — Colonies flat, spreading, cover- 
ing dish after 2 wk, with moderate aerial mycelium, and even, 
lobate margins. On MEA surface luteous buff, reverse orange; 
on PDA surface luteous buff, reverse saffron; on OA surface 
luteous buff. 


Notes — Li et al. (2015) described several species of Para- 
stagonospora from Dactylis collected in Italy, namely P. allo- 
uniseptata (conidia 1-septate, 16—22 x 2.5-3.5 um), P. dac- 
tylidigena (conidia 3-septate, 25—40 x 4-5.5 um), P. minima 
(conidia 3-septate, 20—28 x 3.5—4.5 um) and D italica (conidia 
3-septate, 25-32 x 3—4 um). An additional three sexual spe- 
cies were subsequently described from this host, namely 
P. campignensis (Li et al. 2016), P. fusiformis and P. poaceicola 
(Thambugala et al. 2017). Parastagonospora dactylidigena 
which also occurs on Dactylis, is phylogenetically distinct, 
and morphologically easily distinguished from these taxa and 
those in the P. nodorum complex, based its multiseptate, longer 
conidia. 


Parastagonospora golestanensis B.A. McDonald, P.C. 
Brunner, Croll, D. Pereira & Crous, sp. nov. — MycoBank 
MB838013; Fig. 7 


Etymology. Name reflects the location where it was collected, Golestan 
Province, Iran. 


Typus. IRAN, Golestan Province, on Agropyron tauri, 2011, M. Razavi (holo- 
type CBS H-24388, cultures ex-type CPC 36217 = Pat-6 IR_6_1.1 = CBS 
146871). 


Conidiomata solitary, pycnidial, brown, globose, 200-350 um 
diam, with central ostiole oozing saffron conidial masses; wall 
with 3—6 layers of pale brown textura angularis. Conidiophores 
reduced to conidiogenous cells lining the inner cavity. Conidio- 
genous cells hyaline, smooth, ampulliform, 5—10 x A-B um, 
proliferating percurrently at apex. Conidia solitary, hyaline, 
smooth, guttulate, subcylindrical, straight to slightly curved, 
apex subobtuse, base truncate, 1.5—2 um diam, (1—)3-septate, 
(22-)25-30(-35) x 2.5(—3) um. 


Fig.6 Parastagonospora dactylidigena (CPC 36213). a. Conidiomata developing on PDA; b—c. conidiogenous cells giving rise to conidia; d. conidia. — Scale 


bars: a—c = 200 um, all others = 10 um. 


Fig. 7 Parastagonospora golestanensis (CPC 36217). a. Conidiomata developing on PDA; b—c. conidiogenous cells giving rise to conidia; d. conidia. — Scale 


bars: a = 350 um, all others = 10 um. 
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Culture characteristics — Colonies flat, spreading, reaching 
60 mm diam after 2 wk, with abundant aerial mycelium, and 
even, lobate margins. On MEA surface saffron, reverse sienna; 
on PDA surface and reverse saffron; on OA surface saffron. 


Notes — Parastagonospora golestanensis occurs on Agro- 
pyron, but has also been isolated from Dactylis glomerata. 
Morphologically, it is similar to taxa occurring in the P. nodorum 
complex, and these species are best distinguished based on 
DNA data. 


Parastagonospora jasniorum B.A. McDonald, P.C. Brunner, 
Croll, D. Pereira & Crous, sp. nov. — MycoBank MB838014; 
Fig. 8 


Etymology. Name reflects a condensation of Jasmine and Nicholas, the 
first names of the children of the senior author Patrick Brunner. 


Typus. Iran, Golestan Province, near Aliabad-e-Katul, on infected glume 
of Triticum aestivum, June 2005, M. Razavi (holotype CBS H-24383, cultures 
ex-type CPC 36200 = P-2 IR_A1_3.1A= CBS 146866). 


Ascocarps pseudothecial, solitary to aggregated, 120—170 um 
diam, globose, brown, with central ostiole; wall of 3—6 layers of 
brown textura angularis. Asci 55—70 x 8—11 um, bitunicate, short 
stipitate, straight to flexuous, subcylindrical, with well-defined 
apical chamber, 1—1.5 um diam. Pseudoparaphyses hyaline, 
septate, hyphae-like, branched, septate, 2.5-3 um diam. 
Ascospores bi- to tri-seriate, fusoid, pale brown, verruculose, 
widest in penultimate cell, 3-septate, not to slightly constricted 
at septa, (19—)21-22(—23) x (3.5-)4(—5) um. Conidiomata 
250-300 um diam, dark brown to black, pycnidial, globose 
with central ostiole, exuding pale brown conidial mass; wall of 
2-3 layers of brown textura angularis. Conidiophores reduced 
to conidiogenous cells lining the inner cavity. Conidiogenous 
cells phialidic, hyaline, smooth, aggregated, ampulliform, with 
percurrent proliferation at apex, 5—6 x 4-5 um. Conidia pale 
brown, smooth, straight to flexuous, subcylindrical, guttulate, 
apical cell with slight taper to subobtuse apex, basal cell with 


slight taper to truncate hilum, 1.5-2 um diam, (1—)3(—5)-sep- 
tate, (22—)27—32(-35) x (2.5—)3 um. 

Cultural characteristics — Colonies covering dish in 2 wk at 
25 °C in the dark, with moderate fluffy aerial mycelium and 
smooth, lobate margin. On MEA surface hazel with isabelline 
outer region, reverse cinnamon with patches of isabelline; on 
PDA surface and reverse honey; on OA surface honey with 
patches of hazel. 


Notes — Parastagonospora jasniorum is distinguished from 
P. nodorum (conidia 0—3-septate) based on its conidial septa- 
tion. 


Parastagonospora nodorum (Berk.) Quaedvlieg et al., Stud. 
Mycol. 75: 363. 2013 — Fig. 9 


Basionym. Depazea nodorum Berk., Gard. Chron., London: 601. 1845. 

Synonyms. Septoria nodorum (Berk.) Berk., Gard. Chron., London: 601. 
1845. 

Stagonospora nodorum (Berk.) E. Castell. & Germano, Annali Fac. Sci. 
Agr. Univ. Torino 10: 71. 1977. [1975—1976]. 

Leptosphaeria nodorum E. Mull., Phytopath. J. 19: 409. 1952. 

Phaeosphaeria nodorum (E. Mull.) Hedjar., Sydowia 22: 79. 1969 [1968]. 


Conidiomata solitary, pycnidial, brown, globose, 200-300 um 
diam, with central ostiole, 10-15 um diam, oozing hyaline 
conidial masses; wall with 3—6 layers of pale brown textura 
angularis. Conidiophores reduced to conidiogenous cells lining 
the inner cavity. Conidiogenous cells hyaline, smooth, globose 
to ampulliform, 5-7 x 4—6 um, proliferating percurrently at 
apex. Conidia solitary, hyaline, smooth, guttulate, subcylindrical, 
straight to irregularly curved, apex subobtuse, base truncate, 
1—3-septate, (11-)13—18(—28) x (2.5—)3(—3.5) um. 

Culture characteristics — Colonies flat, spreading, cover- 
ing dish after 2 wk, with moderate aerial mycelium, and even, 
smooth margins. On MEA, PDA and OA surface pale olivaceous 
grey to olivaceous grey, reverse olivaceous grey. 
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Fig. 8 Parastagonospora jasniorum (CPC 36200). a. Conidiomata developing on PDA; b—c. conidiogenous cells giving rise to conidia; d conidia; e—h. asci 


and ascospores. — Scale bars: a = 300 um, all others = 10 um. 
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Material examined. IRAN, Golestan Province, near Aliabad-e-Katul, on Triti- 
cum aestivum, June 2005, M. Razavi (neotype specimen designated here CBS 
H-24391, MBT394790, culture ex-neotype CPC 36202 = nodorum_IR_B2.1B 
= CBS 146873). 


Notes — Bissett (1982) noted that conidia in vivo can be 

1—3-septate, 13—28 x 2.8—4.6 (av. 19 x 3.6) um in size, while 
Eyal et al. (1987) reported conidia to be 0O—3-septate, 15—32 x 
2—4 um. 
Parastagonospora nodorum is a common pathogen of Triticum 
aestivum, and occurs wherever this host is grown. This patho- 
gen was originally described as Depazea nodorum from wheat 
collected in the UK, but no holotype specimen was indicated. 
This is rectified here, with the present reference strain (for 
which a full genome sequence is available) being designated 
as ex-neotype. 


Parastagonospora pseudonodorum B.A. McDonald, P.C. 
Brunner, Croll, D. Pereira & Crous, sp. nov. — MycoBank 
MB838015; Fig. 10 


Etymology. Name reflects its morphological similarity to Parastagonospora 
nodorum. 


Typus. Iran, Golestan Province, on infected Triticum aestivum, 2010, 
M. Razavi (holotype CBS H-24384, cultures ex-type CPC 36208 = Pat-1 
IR_5.2B = CBS 146867). 


Conidiomata solitary, pycnidial, brown, globose, 200—350 um 
diam, with central ostiole; wall with 3—6 layers of pale brown 
textura angularis. Conidiophores reduced to conidiogenous cells 
lining the inner cavity. Conidiogenous cells hyaline, smooth, 
ampulliform to subcylindrical, 4—9 x 4—6 um, proliferating per- 
currently at apex. Conidia solitary, hyaline, smooth, guttulate, 
cylindrical, straight to slightly curved, apex subobtuse, base 
truncate, 2.5-3 um diam, 3-septate, (27—)30-—33(—36) x (2.5-) 
3(-3.5) um. 


Culture characteristics — Colonies spreading, covering dish 
after 2 wk, with moderate fluffy aerial mycelium, and even, 
lobate margins. On MEA surface saffron, reverse orange; on 
PDA surface and reverse saffron; on OA surface smoke grey. 


Notes — Parastagonospora pseudonodorum coexists with 
P. nodorum on wheat, but is more often found on heads and 
seeds instead of on leaves. Morphologically the two species 
can be distinguished in that conidia of P pseudonodorum are 
somewhat longer than those of P. nodorum (see above). 


Parastagonospora stipae B.A. McDonald, P.C. Brunner, 
Croll, D. Pereira, Mordecai & Crous, sp. nov. — MycoBank 
MB838016; Fig. 11 


Etymology. Name reflects the host genus it was isolated from, Stipa. 


Typus. USA, California, Woodside, Jasper Ridge Biological Preserve, on 
Stipa pulchra (formerly Nassella pulchra), 2015, E. Mordecai and E Spear (holo- 
type CBS H-24389, cultures ex-type CPC 36223 = New-3 MORDECAI_ 1418 = 
CBS 146872). 


Culture nearly sterile, with only a few conidiomata observed. 
Conidiomata solitary, pycnidial, dark brown, globose, 150— 
180 um diam, with central, darker brown ostiole; wall with 3—6 
layers of pale brown textura angularis. Conidiophores reduced 
to conidiogenous cells lining the inner cavity. Conidiogenous 
cells hyaline, smooth, ampulliform to subcylindrical, 5-6 x 
3-4 um, proliferating percurrently at apex. Conidia solitary, 
hyaline, smooth, guttulate, subcylindrical, straight to slightly 
curved, apex subobtuse, base truncate, 2 um diam, 1-septate, 
(8—)10—13(-18) x (2.5—)3 um. 

Culture characteristics — Colonies flat, spreading, cover- 
ing dish after 2 wk, with moderate aerial mycelium, and even, 
smooth margins. On MEA surface luteous, reverse sienna; on 
PDA surface and reverse sienna; on OA surface pale luteous. 


Fig. 9 Parastagonospora nodorum (CPC 36202). a. Conidiomata developing on PDA; b. conidial cirrhus; c. conidiogenous cells giving rise to conidia; 


d. conidia. — Scale bars: a = 300 um, all others = 10 um. 


Fig. 10 Parastagonospora pseudonodorum (CPC 36208). a. Conidiomata developing on PDA; b—c. conidiogenous cells giving rise to conidia; d. conidia. — 


Scale bars: a = 350 um, all others = 10 um. 
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Notes — Parastagonospora stipae occurs on Stipa, and is 
distinguished from taxa in the P. nodorum complex based on 
its smaller conidia. 


Parastagonospora zildae B.A. McDonald, P.C. Brunner, Croll, 
D. Pereira & Crous, sp. nov. — MycoBank MB838017; Fig. 
12 


Etymology. zildae = Named after Zilda, the first name of the mother of 
the co-author Danilo Pereira. 


Typus. Iran, Golestan Province, near Aliabad-e-Katul, on infected glume 
of Triticum aestivum, June 2005, M. Razavi (holotype CBS H-24381, cultures 
ex-type CPC 36198 = P-1 IR_B4.2A = CBS 146864). 


CPC 36198 = P-1 IR_B4.2A: Ascocarps pseudothecial, solitary 
to aggregated, erumpent, globose, dark brown, 200—250 um 
diam; wall of 3—6 layers of brown textura angularis; ostiole 
central, 15—25 um diam, not to slightly papillate. Asci 55-85 
x 7—9 um, bitunicate, subcylindrical, straight to slightly curved, 
short stipitate, 8-spored, with well-defined apical chamber, 
1-1.5 um diam. Pseudoparaphyses filiform, hyaline, septate, 
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hyphae-like, unbranched, rarely anastomosing, 2-2.5 um 
diam. Ascospores bi- to triseriate in asci, fusoid, guttulate, pale 
brown, smooth, widest in penultimate cell, becoming slightly 
constricted at septa with age, (24—)26—28(—32) x (3.5—)4 um, 
5-6-septate. 
CPC 36221 = New-1 IR_C2.2A: Ascocarps pseudothecial, 
solitary, erumpent, globose, dark brown, 180—250 um diam; 
wall of 3—6 layers of brown textura angularis; ostiole central, 
15-25 um diam, not to slightly papillate. Asci bitunicate, 
65-85 x 8—11 um, subcylindrical, straight to slightly curved, 
short stipitate, 8-spored, with well-defined apical chamber, 
1.5-2 um diam. Pseudoparaphyses filiform, hyaline, septate, 
hyphae-like, branched below, rarely anastomosing, 2—3 um 
diam. Ascospores bi- to triseriate in asci, fusoid, guttulate, pale 
brown, smooth, widest in penultimate cell, becoming slightly 
constricted at septa with age, (24—)25—26(—28) x 4(—5) um, 
5(—6)-septate. 

Cultural characteristics — Colonies covering dish after 2 wk 
at 25 °C in the dark, smooth, with even margin, and moderate 
to fluffy aerial mycelium. On MEA surface buff with patches 


2 2f SS e $ T t 
Fig. 11 Parastagonospora stipae (CPC 36223). a. Broken conidiomata with conidia; b—c. conidiogenous cells giving rise to conidia; d. conidia. — Scale bars: 
a = 200 um, all others = 10 um. 
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Fig. 12 Parastagonospora zildae (CPC 36198). a-c. Ascomata developing on PDA; dh asci, pseudoparaphyses and ascospores; i. ascospores. — Scale 


bars: a—c = 200 um, all others = 10 um. 
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of honey, reverse isabelline; on PDA surface pale olivaceous 
grey, reverse grey olivaceous; on OA surface grey olivaceous 
in centre, cinnamon in outer region. 


Additional specimen examined. IRAN, Golestan Province, near Aliabad-e- 
Katul, on infected glume of Triticum aestivum, June 2005, M. Razavi (CBS 
H-24382, cultures CPC 36221 = New-1 IR_C2.2A = CBS 146865). 


Notes — Cultures of P. zildae only produced a sexual morph 
in culture, making a morphological comparison with asexual 
species difficult. Species such as P. nodorum do have a sexual 
morph (originally described as Leptosphaeria nodorum, asco- 
spores 3-septate, 23-32 x 4—6 um; Eyal et al. 1987), although 
this is not commonly observed. 


DISCUSSION 


Genome-scale sequence comparisons allowed us to differenti- 
ate several new species of Parastagonospora and determine 
their phylogenetic relationships. The highest species diversity 
was found in Parastagonospora collections from the Fertile 
Crescent, the region where wheat was domesticated, consistent 
with earlier hypotheses that the Fertile Crescent is the centre of 
origin for wheat-infecting Parastagonospora spp. We consider 
it noteworthy that five Parastagonospora species were found 
infecting wheat in Iran, but only two of these species (P. no- 
dorum and P. pseudonodorum) are commonly found on wheat 
worldwide. This finding suggests that there is potential for other 
Parastagonospora species to ‘escape’ from the Fertile Crescent 
region and emerge as new pathogens on wheat in other regions. 
Both P. nodorum and P. pseudonodorum are known to infect 
seed, so it is likely that these pathogens were moved around 
the world on infected seed, but it remains unknown whether 
the other three Parastagonospora species found on wheat in 
Iran are also able to infect seed. 


An earlier coalescent analysis of five genes suggested that the 
grass-infecting Parastagonospora species could be broadly 
separated into ‘young’ and ‘old’ species, with the youngest spe- 
cies including P. pseudonodorum, P. golestanensis, P. nodorum 
and D avenae and the older species including P. zildae and 
P. jasniorum (McDonald et al. 2012). The new genome-scale 
analyses largely supported this hypothesis, with the younger 
species showing shorter branch lengths consistent with more 
recent speciation events while the older species showed deeper 
divergences consistent with much earlier speciation events. We 
speculate that the young Parastagonospora species emerged 
mainly in the Fertile Crescent as a result of selection to become 
specialized to infect agricultural crops (e.g., wheat, oats and 
barley) and the weedy grasses often associated with these 
crops. The most parsimonious scenario is that the speciation 
events coincide with the onset of agriculture in the region, 
similar to findings on the emergence of the wheat pathogen 
Zymoseptoria tritici (syn. Mycosphaerella graminicola) (Stuken- 
brock et al. 2007a). 


An earlier study (McDonald et al. 2012) that analyzed se- 
quences from five genes in more than 350 global Parastagono- 
spora strains identified 14 strains that appeared to be hybrids 
because they carried alleles that were otherwise exclusive to 
either P. nodorum or P. pseudonodorum (formerly Pat-1). The 
complete genome sequences reported here for three of these 
strains showed no evidence for hybridization, with each isolate 
falling clearly into either the P nodorum (CASSILS) or the 
P. pseudonodorum (JANSEN4, HARTNEY) clade. Genomes 
of true hybrids are expected to occupy an intermediate phylo- 
genetic position compared to the parental species. Future 
analyses of the genomic datasets should focus on identifying 
rare recombination events that can lead to gene introgres- 
sion among species. We expect that such events could have 


occurred among P pseudonodorum and P. nodorum, which 
could also explain how P. pseudonodorum acquired the ToxA, 
Tox1 and Tox3 genes from P. nodorum (McDonald et al. 2013, 
Ghaderi et al. 2020). The preponderance of P. pseudonodorum 
on grain and ears suggests that this species exhibits host-tissue 
specialization, with a preference to infect ears. In comparison, 
P. nodorum is typically found on both leaves and ears. 


Overall, our data is consistent with the hypothesis that P. no- 
dorum originated as a pathogen of wild grasses in the Fertile 
Crescent, then emerged as a wheat pathogen via host-tracking 
during the domestication of wheat in the same region. It is likely 
that P. nodorum became distributed globally via movement of 
infected wheat seed during the global spread of wheat agricul- 
ture that included North America. The discovery in California of 
Stipa pathogens that appeared closely related to P. nodorum 
is also consistent with a ‘spill-back’ as a possible source for 
this population infecting wild grasses. While the tree based on 
five AFTOL genes provides some weak support for the spill- 
back hypothesis, the tree based on four AFTOL genes does 
not. Therefore, we do not find definitive evidence that P. stipae 
was most closely associated with P. nodorum or consistently 
shared a most recent common ancestor with P. nodorum. Ge- 
nome sequences from P. stipae will be needed to further test 
the spillback hypothesis. 
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Appendix Maximum likelihood phylogenetic tree based on ITS sequences of 26 known Parastagonospora species. Additional ITS sequences were obtained 
from Goonasekara et al. (2019), Marin-Felix et al. (2019) and Brahmanage et al. (2020). 


